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1. Introduction 

In lattice simulations of QCD on a cubic volume (V = L 3 ) with periodic boundary conditions 
imposed on the fields, hadronic momenta, p, are quantized according to pi = 2n/L x n,-, where 
i = 1,2,3 and the n\ are integers. This implies that on currently available lattices the gaps between 
neighbouring momenta are large, typically about 500 MeV or so, limiting the phenomenological 
reach of the simulations. In ref. [[I], 1 Bedaque proposed the use of twisted boundary conditions 
on the quark fields (q(x)), e.g. 

q(xi+L) = e Wi q(xi), (1.1) 

so that the momentum spectrum is 

Thus by varying the twisting angles, the 0,'s, arbitrary momenta can be reached. In refs. [^, ^ we 
have investigated the use of twisted boundary conditions theoretically and numerically and in this 
note we briefly report the conclusions of these studies. 



2. Finite- Volume Effects with Twisted Boundary Conditions 

The three main results from ref. [[|] are: 

1. For physical quantities without final state interactions, such as masses or matrix elements 
of local operators between states consisting of the vacuum or a single hadron, the flavour 
symmetry breaking induced by the twist only affects the finite-volume corrections, which 
nevertheless remain exponentially small. 

2. For physical quantities without final state interactions, one can also use partially twisted 
boundary conditions in which the sea quarks satisfy periodic boundary conditions, but some 
or all of the valence quarks satisfy twisted boundary conditions, with exponential precision 
in the volume. This implies that in unquenched simulations it is not necessary to generate 
new gluon configurations for every choice of twisting angle, making the method much more 
practicable. 

3. For amplitudes which do involve final state interactions (such as K —> %% decays), it is not 
possible in general to extract the physical matrix elements using twisted boundary conditions 
(at least without introducing new ideas) 

We now briefly discuss each of these results in turn. 

2.1 Result 1. 

The choice of boundary conditions for the fields only affects the finite-volume corrections 
and when extracting physical quantities we need to establish that these corrections are negligibly 
small. Infrared effects dominated by the pseudo-Goldstone bosons can be evaluated using Chiral 
Perturbation Theory and in ref. [||] we derive the chiral Lagrangian corresponding to the twisted 
boundary conditions. Consider SU (N)l x SU (N)r chiral symmetry and impose the twisted bound- 
ary conditions on the quark fields q{xj +L) = Uiqixi) = exp(idf , T a )q(xi), where we take the £/; 

1 See also the references cited in [Nl 0] for earlier related ideas. 
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to be diagonal N x N unitary matrices and the T"'s are the generators of the Cartan subalgebra 
of the (light) flavour U (N)v group. In the chiral Lagrangian these boundary conditions imply 
that T(x), the coset representative of SU (N)l x SU (N)r/SU (N)y satisfies the boundary conditions 
E(x, +L) = UiL(xt)uJ . The Feynman rules for the Chiral Lagrangian are the usual ones, except 
that the momentum spectrum for the Pseudo-Goldstone bosons is that expected of a particle com- 



posed of a quark-antiquark pair, each with a momentum given in terms of its twist by eq. (|1.2|). The 
effect of the twist is a change in the charged meson's dispersion relation, 



J Tt,p 



m it,p + 



/ 27T7? - d 



(2.1) 



where m n j p is the meson mass and 0, = 8 u j — dd.i, the difference of the twists of the u and d quarks. 

As an example consider the finite- volume corrections in f n , the leptonic decay constant of the 
7T-mesons. As L — > oo the leading corrections are found to be: 



-m„L 



fl (2nm K L)V2 \ 3 



-£cOS0 + l 



i=l 



-m„L 



fl (2nm K L)yi \ 3 f 



^£cos0,- ] (2.2) 



where A(X)/X = (X(L) — X(oo))/X(oo). The results in eq.( |2.2[ ) illustrate the point that whereas 
the finite-volume corrections do depend on the boundary conditions as expected, they nevertheless 
remain exponentially small (see ref. ^ for further examples). The difference of the correction 



terms in eq.(2.2) is a manifestation of the violation of isospin by the boundary conditions. 



2.2 Result 2. 

The derivation of the Chiral Lagrangian and the corresponding Feynman rules with Partially 
Twisted boundary conditions follows similar steps to that for partially quenched QCD. As an illus- 
tration of the results consider the finite- volume corrections to fy, the leptonic decay constant of the 
,K-meson, with the d and s quarks satisfying periodic boundary conditions and the u quark with (a) 
untwisted, (b) fully twisted (c) partially twisted boundary conditions: 



fx* 



"4 71" (2^1)3/2 



1 n I?-, COS 0,-f) 



ml e - m * L (y3 pn „ a 3\ 



(a) 
(b) 
(c) 



(2.3) 



This example illustrates again that, in general, the finite-volume corrections are different for the 
three cases but they are always exponentially small. 



2.3 Result 3. 

To illustrate the problem in using twisted boundary conditions for processes with final state 
interactions consider K —> Tin decays in the 1 = channel. The twisted boundary conditions break 
isospin symmetry which implies that energy eigenstates are no longer states with definite isospin. 



352/3 



Twisted Boundary Conditions in Lattice Simulations 



Christopher Sachrajda 



This is manifest even in the free theory. As an example consider the two pions to be in the centre-of- 
mass frame and let the u and d quark satisfy twisted and periodic boundary conditions respectively 
so that the lowest momenta of the neutral and charged pions are p^o = and p n ± = ±d/L. The 
corresponding energies of the two-pion states are then E^o^o = 2m K and E w + % - = 2^m\ + d 2 /L 2 
making it explicit that the energy and isospin eigenstates are different. 

The interacting theory contains 7T + 7r~ <-> % % transitions which complicate the analysis very 
significantly making it impossible to relate the finite-volume energy eigenstates to the infinite- 
volume energy and isospin eigenstates \%%) i= q and \%1t)i-2- Without new ideas, it not therefore 
not possible to determine physical / = K —> %% amplitudes using twisted boundary conditions 
with finite-volume corrections under control. 

There is a further difficulty which we would like to exhibit. Even if we were able to overcome 
or circumvent the problem above (for example, by using G-parity boundary conditions proposed 
by Christ and Kim [||] 2 ). The finite-volume effects which decrease only slowly, i.e. as powers of 
the volume, come from the propagation of two-pion states illustrated by the diagram: 

P 



The circles represent insertions which are two-particle irreducible in the s-channel. In evaluating 
finite-volume corrections we replace the sums over the discrete momenta in a finite volume by 
integrals over continuous infinite-volume momenta using the Poisson summation formula: 

d 3 p 



1 



E/(p) 



I 

fez 3 



(2k) 



(2.4) 



For two-particles propagating with an energy above the two-pion threshold (E > 2m n ) there are 
poles in f(p) (corresponding to the cut in infinite volume) and this implies that the terms with 
I ^ in eq. (2A) are not exponentially small. For I = K —> Tin decays, the sea quarks contribute 
to the final-state interactions as illustrated for example, by the following diagram: 




Therefore, unless the same boundary conditions are imposed for both the valence and the sea 
quarks, there is a breakdown of unitarity induced by the partially twisted boundary conditions 
and the finite-volume corrections are not under control. This is analogous to the inconsistencies 
encountered in quenched and partially quenched QCD for K ^ %% decays above the two-pion 
threshold in refs. [|6[ \j§ . 



3. Numerical study of partially twisted boundary conditions 

This section describes the results of an exploratory numerical study of partially twisted bound- 
ary conditions, in which we investigate, whether the meson's momentum follows eq.(2T) and 



"The twist angle is now fixed to be n and so the benefit of using partial twisting are not so clear. 
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whether the meson decay constants f K and f p and the matrix element Zp = (0\P\n) of the pseudo 
scalar density P are independent of the twisting angle. 



The study has been carried out on 200 independent, non-perturbatively improved Nf = 2 Wil- 
son gauge configurations with the following specifications: j8 = 5.2, (L/a) 3 xT/a = 16 3 x 32, 

a ~ 0. lfm, JCyal = K"sea = 

0.13500 and 0.13550 corresponding to m„/m p =0.697(11) and 0.566(16) 

(see also 

We evaluated mesonic correlation functions (details in ref. [Q]) for all possible pairs of valence 
quark twisting angles 6 uc i G {0, (2,0,0), (0, 71,0) and (3,3,3)} at the Fourier momenta p\ at = and 
(0,±27r/L,0). In addition we also evaluated all quantities at \p\ at \ = V2 x 2n/L and v3 x 2n/L 
for vanishing twist. The jack-knife analysis of the resulting correlation functions focused on the 
momentum dependence of E n and E p and we also monitored the un-improved bare decay constants 
f n and f p and matrix element Zp under the variation of the momentum. 

Here we discuss only the data for the set at m n jm p = 0.697(11) and refer the reader to ref. 
[Q] for details and qualitatively equivalent results for the lighter data set. The plots in figure |] 
show our results as a function of (pL) 2 . The positions of the discrete Fourier momenta \p\atL\ = 0, 
2n, y/2 x 2n and y/3 x 2n which can be reached without twisting are indicated by dashed vertical 
lines. In each plot the (blue) triangles correspond to points in which the correlation function was 
evaluated with p\ at = 0, but with all possible pairs of d Ll and 8d- The (red) diamonds and (green) 
squares represent the results obtained with p\ &x = (0,27r/L,0) and p\ at = — (0,2n/L,0) respectively, 
combined with all possible pairs of 6 Ll and 6^. The four points with 8 U = 8d = but \pm\ = 
0, 2n/L, y/2 x 2n/L and y/3 x 2n/L are denoted by (black) circles. The three main results are: 



1. The energies of n and p -mesons (with masses below the two-pion threshold) satisfy eq. ([2JJ) 
very well, particularly for small values of the momentum where lattice artifacts are small. 

2. The values of the leptonic decay constants f K and f p and of the matrix element Zp are in- 
dependent of the twisting angles for small values of the momentum. Deviations for large 
momenta are of the same magnitude when using twisted boundary conditions or only Fourier 
transformation to induce momentum. 

3. Twisted boundary conditions do not introduce additional noise in the data. The combined 
statistical and systematic error on the meson masses and matrix elements increases smoothly 
when increasing the meson's momentum by varying the angles 6 U ^- However, when com- 
paring results obtained with twisted and periodic boundary conditions with similar momenta 
(i.e. around 2n/L or y/2(2n/L)), the errors are found to be comparable. 

4. Conclusions 



The theoretical and numerical results presented in this note demonstrate that the use of partially 
twisted boundary conditions does indeed improve the momentum resolution in lattice phenomenol- 
ogy at relatively little cost for physical quantities without final-state interactions. We now look for- 
ward to applying these techniques to studies of weak matrix elements and hadronic structure. We 
have also demonstr ated that twisted and partially twisted boundary conditions cannot be applied in 
general to physical quantities with final state interactions such as K — > nn decays with isospin 0. 
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m n /m p = 0.70 



0.12r 



m n /m p = 0.70 




Figure 1: Results for the case m n /m p = 0.70. The plots in the first line illustrate the results for the dispersion 
relation for the % and the p (empty and full symbols respectively) and the associated error in the case of the 
n as a function of the momentum. The second line shows the results for the % and p decay constant and Zp. 
In each plot the horizontal lines represent the central value at p\ A = Q\ = 62 = 0. Note also in all the four 
plots the superimposed data points at |pi at | = 0, 2n/L, which correspond to the cases Q„ — C / = 0. 
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